Construction of model dielectric functions for two and three dimensional electron 

liquids from density functionals 
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The Thomas-Fermi (TF) approximation for the static dielectric function of a three-dimensional 
(3d) electron liquid can be derived by minimizing the TF local-density approximation for the kinetic- 
energy functional. Here we show that this connection between energy functionals and model dielec- 
tric functions is not an artifact, but a general paradigm. Four examples are worked out in detail, 
by calculating the model dielectric functions that follow, respectively, from (i) exchange corrections 
to TF theory in 3d, i.e., TF-Dirac theory, (ii) further correlation corrections to TF-Dirac theory in 
3d, (iii) TF theory in 2d, and (iv) exchange corrections to TF theory in 2d. Each of these cases has 
certain interesting features, revealing connections between independent many-body methods, that 
are discussed in some detail. As a byproduct of these investigations we also find that a common 
textbook statement about the long-wavelength (k — * 0) limit of the random-phase approximation is 
not fully correct. 
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I. INTRODUCTION 

Screening is one of the most important manifesta- 
tions of many-body effects in metals, semiconductors and 
plasmas. 1,2 In the Thomas-Fermi (TF) approximation, 3 
static screening in three-dimensional (3-d) electron liq- 
uids is described by means of a wave-vector-dependent 
dielectric function, 4-6 



e TF (k) = 1 



k 2 

_ fc 2_: 



(1) 



where the TF screening wave vector is krF = y4fep /ircto- 
Here fcp = (Stt 2 ^) 1 ^ 3 and ao = ft 2 /me 2 denote the Fermi 
wave vector and the Bohr radius, respectively, and no is 
the charge density of the unperturbed electron liquid. 
This model dielectric function is sometimes used in its 
own right as an approximation for screening in metals, 
plasmas and doped semiconductors; more frequently, it 
appears as starting point for, or limiting case of, more 
sophisticated many-body treatments of screening. 4-8 We 
stress from the outset that our aim in the present work is 
not the construction of high-precision or material-specific 
expressions for the dielectric function, but to explore con- 
nections between different many-body methods that be- 
come visible through the expressions they yield for the 
dielectric functions. 

The Thomas-Fermi approximation to the kinetic en- 
ergy of a 3-d electron liquid is 9 ' 10 



3h 2 (3ir 2 ) 2 / 3 



10m 



An(r) 5/3 = C / d 3 m(r) 5/3 



T TF [n] = 

Willi I I 

(2) 

While both exF (k) and 7V^[n] are widely known and 
employed concepts, it is worthwhile stressing that the 
two uses of the label 'Thomas-Fermi' are not equivalent: 



the 3-d TF dielectric function can be obtained from the 
TF energy functional by one additional approximation, 
a linearization of the Eulcr equation resulting from min- 
imizing T TF [n]. 7 ' n To make our paper self-contained we 
start, in Sec. II, by briefly recalling this derivation. 

The inclusion of many-body effects beyond the TF 
approximation is in general highly nontrivial, and typ- 
ically accomplished by either constructing local-field cor- 
rections to the random-phase approximation (RPA), 7 or 
by first-principles density-functional calculations of the 
dielectric functions of realistic materials. 12 In Sec. Ill we 
show that certain important many-body effects can be 
incorporated into the TF dielectric function in a much 
simpler way, by employing the Dirac-Slater local-density 
approximation (LDA) for the exchange energy, 8-10 



Ei DA \n] 



3e^ / 3 
4 Itt 



1/3 



d 3 r m 



(r )4/3 

D J d 3 rn{vf' 3 . 



(3) 
(4) 



Addition of this to the TF energy functional 
leads to the so-called Thomas-Fermi-Dirac (TFD) 
approximation. 8-10 Further correlation corrections can 
be described by the local-density approximation 9 



EL DA [n] 



Ex DA [n]+ I d 3 re c {n)\ n ^ n{r) , (5) 



where e c (n) is the per- volume correlation energy 
of the uniform electron liquid 13 underlying common 
parametrizations of the LDA. 9 In Sec. V we derive the 
dielectric function arising from £^ c . Sees. VI and VII 
extend these calculations to two-dimensional electron liq- 
uids, within TF and TFD theory, respectively. As a 
byproduct of this investigation we also revisit, in Sec. IV, 



2 



the long wavelength limit of the three-dimensional dielec- 
tric function found in the RPA, which has as nontrivial 
correction similar to the one found here within TFD the- 
ory. 



II. TF DIELECTRIC FUNCTION IN 
THREE-DIMENSIONAL ELECTRON LIQUIDS 

This derivation, although elementary, will be presented 
in some detail, because it serves as a model for the more 
complicated cases treated afterwards. The 3-d TF total- 
energy functional is 



E TF [n]=T TF [n] 

?/" 3r / dV T^ + / rf3rn(r)v - t(] 



+ 

*(r), (6) 



where the second term on the right-hand-side (rhs) is the 
electrostatic (Hartree) energy of the charge distribution 
n(r) with itself, and the third term describes the poten- 
tial energy of this charge distribution in the external po- 
tential v ex t (r) , composed of a smeared-out background of 
positive charge with constant density nb g and potential 
Vb g and an additional test charge n tes t( r ) with potential 
^test( r )- It is the screening of this test charge, which is 
normally taken to be a point charge, that the theory aims 
to describe. Equilibrium is characterized by the condi- 
tion S(ETF[n] — nN)/5n(r) = 0, where the total particle 
number N = j d 3 rn(r) and the chemical potential u at 
zero temperature is the Fermi energy tp. This minimiza- 
tion leads to the Euler equation 



Cn(r) 2 / 3 + v ext (r) + v H {r)-fi = 0, 



(7) 



where the Hartree potential is vh = SEH[n]/Sn = 
e 2 J dVra(r')/|r - r'\. Solving Eq. (7) for n(r) we find 



3/2 



(8) 



where we have introduced the effective, or screened, po- 
tential v s = vh + v ext . If this potential is weak, i.e., 
v s /^i <§C 1, Eq. (8) can be linearized to give 



n(r) 



3/A 
5C7 



3/2 



3 v s (r) 
2 ii 



(9) 



By introducing the charge density of the unperturbed 
system, n — k F /3n 2 (where ft 2 hp /2m = fj,), and the 
corresponding 3-d density of states, g a — mk F /h 2 n 2 , we 
can write this as n(r) = uq — gov s (r) =: uq + rii n d- The 
screened potential becomes 



v s {r) = v ext (r)+e 2 J d 3 r 



3 , n (r') + n ind (r') 



(10) 



Stability of the unperturbed electron liquid (with v tes t = 
0) requires no = —ni, g , so that the background- 
contribution to v ex t is cancelled precisely by the first term 
arising from the integral. Hence, 



V s (r) = Vtest 



(r) + e 2 J d' 



.3 / n ind( r ') 



Upon Fourier transformation this becomes 

v s (k) = v t est(k) + e 2 J r [n m d(r)]J c "[l/r], 



(11) 



(12) 



where T stands for the 3-d Fourier transform opera- 
tor. In the second term we have used the convolu- 
tion theorem to decompose the Fourier transform of the 
Hartree potential into a product of two Fourier trans- 
forms. By substituting rii n d = —gov s in this equation, 
using F[l/r] — 4ir/k 2 , and solving for v s , we obtain 



Vs(k) 



Vtest(k) 



Vtest (k) _ Vtest (k) 



1 + e 2 g T[l/r] 1 + 



4k f 



e TF (k) 



, (13) 



in agreement with the definitions following Eq. (1). For 
the screening of a point charge with w test (k) = 47re 2 /fc 2 
this becomes 



w.(k) 



47re 2 



k 2 + k 2 



(14) 



TF 



Minimization of the TF kinetic-energy functional fol- 
lowed by a minimization of the resulting Euler equation 
thus reproduces the TF dielectric function 



e TF (k) = 1 + 



4k F 1 
ciott k 2 



(15) 



In the following sections we explore whether this con- 
nection between energy functionals and model dielectric 
functions remains valid in more complex situations, char- 
acterized by additional exchange and correlation terms in 
the functional, and in two dimensions. 



III. TFD DIELECTRIC FUNCTION IN 
THREE-DIMENSIONAL ELECTRON LIQUIDS 

The fact that dielectric functions can be derived from 
energy functionals opens up the possibility to incorporate 
additional physics into a model dielectric function by us- 
ing a more complete energy functional. To illustrate this 
we now deduce an expression that incorporates exchange 
effects into the TF dielectric function. Starting point is 
the 3-d TFD energy functional 



E T fd [n] = E T f [n] 



E^ DA \n] 



(16) 



where the functionals on the rhs are defined in Eq. (6) 
and (4), respectively. The Euler equation corresponding 
to this functional is 

|cn(r) 2 / 3 + v ext (r) + v H (r) - ^n(r) 1 / 3 - M = 0. (17) 
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In accordance with common practice in density- 
functional theory 9 ' 10 we regard the exchange term on the 
rhs as a contribution to the effective potential, now de- 
fined as 

v s = v ext + v H + v x = v ext +v H - (4D/3)n 1/3 . (18) 

The induced density change is defined in terms of this 
potential as before, via n(r) = n — g v s (r) =: n + 
nmd(r), and the counterpart to Eq. (12) becomes 



v s (k) = v test {\i)+e 2 T[n md {r)]J : 



~nn 1/3 ]- (19) 



To evaluate this expression in closed form we linearize 
again, assuming consistently that rii n d <C n^, i.e., the 
self-consistent modification of the density distribution 
due to the test charge is small compared to the density of 
the original unperturbed system. This allows us to write 
n c/ 3 (l + "ind/3no), and the preceding equation 



1/3 



n 



becomes 



v s (k) = VtestQ*-) + e 2 J r [n i „ (1 (r)]^ r 



1/3! 



4D 



9n, 



2/3 



(20) 



The last term can be handled as before, by using m n d — 
—gov sc . Since n is spatially constant, the Fourier trans- 
form of the second-to-last term contributes only at k = 0. 
Solving the preceding equation for v s we find 



V s (k) 



v test {k) - 4DF[nl /3 }/3 



1 



(21) 



where we have substituted the explicit expressions for 
go, n and D in terms of universal functions and the 
Fermi wave vector. To identify the dielectric function 
we consider again the screening of a point charge, with 
vtest = 47re 2 /fc 2 . After a little algebra we find 

where in the last step we have introduced the usual elec- 
tron gas parameters r s = l/a k F a and a — (4/97T) 1 / 3 . 
Eq. (22) illustrates how exchange-corrections modify the 
TF dielectric function, and also illustrates how by start- 
ing from suitable energy functionals additional physics 
can be built into models for e(k). 

Interestingly, the resulting additive exchange correc- 
tion to the TF dielectric function, — ax s /7r, turns out to 
be the same already known to appear as the multiplica- 
tive Hartree-Fock correction to the sound velocity (s), 
compressibility (k) and spin susceptibility (x) of the 3-d 
electron gas: 1,2 



™(k)- £TF (k) = 1- ^ = 1- 1 = 1- f = ~- 

(23) 




FIG. 1: Three-dimensional Thomas- Fermi (dashed curves) 
and Thomas-Fermi-Dirac (continuous curves) dielectric func- 
tion vs. wave vector k/kF, for r 3 = 6 (upper curves) and 
r 3 — 3 (lower curves). For comparison purposes we also in- 
clude the k — > limit of the static RPA dielectric function 
(dotted curves), as given in Eq. (26). For the density depen- 
dence of the TF and TFD dielectric functions see Fig. 3. 



The behaviour of €tfd (k) for two different densities is 
illustrated in Fig. 1 . Exchange is seen to reduce screening 
(i.e., to lead to a dielectric function closer to e = 1 for all 
k) . This behaviour is easily understood because exchange 
tends to keep particles apart, thereby making screening 
less effective. The limits of validity of Eq. (22) are similar 
to that of ordinary TF theory, restricting the use of (22) 
to high and slowly varying densities (small r s and k — > 0) . 



IV. LONG WAVELENGTH LIMIT OF THE 
RANDOM-PHASE APPROXIMATION 

In the random-phase approximation the static dielec- 
tric function is given by 1,2,4-7 



£RPA(k) = 1 



k 2 
~2k? 



1 



z 2 )ln 



l + x 



l-x 



. (24) 



where x = k/2kp- In textbooks 1 ' 2,5 ' 6,11,14 it is usually 
argued that in the limit x — > the term in square brackets 

tends to 1, because, to linear order in x, In 1 -t^ — > 2x, 

' ' ' 1— X 

so that the RPA dielectric function for long wavelength 
density oscillations becomes €RPA{k — ► 0) = erF(k) and 
thus recovers the TF result. (See, e.g., Eq. (5.65) of 
Rcf. 1, p. 105 of Rcf. 11, p. 140 of Ref. 6, Eq. (5.35) of 
Ref. 2, or Eq. (28.9) of Ref. 14.) 

This statement, however, is not fully correct. More 
precisely, it is not consistent in orders of k. To the order 
in k to which TF theory provides an approximation to 
e(k), i.e., constant plus 1/fc 2 , it is not enough to expand 
the logarithm in Eq. (24) to linear order in x. Instead, 
the expansion must be carried to cubic order, which after 
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multiplication by the prefactor ^-(1— x 2 ) also contributes 
a constant term to e RPA (k). The RPA dielectric function 
thus has the k — > expansion 

e RPA (k -+ 0) = ^fc- 2 + f 1 - i^i) fc° + 0(fc 2 )(25) 

= e TF (k)--^ + 0(fc 2 ).(26) 

In the extreme fc — > limit all constants can be neglected 
compared to 1/fc 2 , and RPA and TF theory do indeed 
predict the same divergence. However, already at the 
level of constant corrections to this limit, RPA and TF 
differ, and the TF dielectric function in the form 1 + 
k\ F I k 2 is never obtained as a long wavelength limit from 
the RPA (except in the extreme case r s = 0). 

Interestingly, the additional term — ar s /3ir needed to 
make the k — > expansion of the static RPA consistent 
in orders of k is, up to a factor 1/3 the same additional 
term we found above as an exchange correction to the 
TF dielectric function. 



V. CORRELATION CORRECTIONS TO THE 
THOMAS-FERMI DIELECTRIC FUNCTION 
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0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

k/k p 

FIG. 2: Two-dimensional Thomas-Fermi (dashed curves) and 
Thomas-Fermi-Dirac (continuous curves) dielectric function 
vs. wave vector k/kF, for r 3 — 4 (upper curves) and r s = 1 
(lower curves) . Inclusion of exchange shifts the dielectric func- 
tion to lower values, the shift being larger for lower densities 
(larger r s ), and for fixed density larger in two than in three 
dimensions. For the density dependence of the 2-d TF and 
TFD dielectric functions see Fig. 3. 



Further many-body effects can be included by adding 
the local-density approximation for correlation, Eq. (5) to 
the TFD energy functional. We introduce the derivatives 



and 



v xc [n](r) = 



/xcN(r,r') = 



SE xc [n] 
6n(r) 

S 2 E xc [n] 
5n(r)8n(r') 



(27) 



(28) 



and note that in the LDA, Eq. (5), the exchange- 
correlation kernel f xc is 



/if>](r,r') = 



d 2 e xc (n) 
dn 2 



|n-m(r)tf(r - r') 



(29) 



=:/Jf>]«5(r-r'), (30) 

where the second equality defines the functional f x ® A [n]. 
The Fourier transform of f xt ^ A is simply related to that 
of the static local-field factor Q (k) via 9 



/if>](k) = -e 2 T 



0(k), 



(31) 



where — e 2 T[l/r] is the Fourier transfrom of the Coulomb 
interaction. 

By going through the same steps as before we find 
for the dielectric function corresponding to correlation- 
corrected TFD theory 



e T FDc(k) = 1 + 



4fc, 



+ flb/xc [«o 



(32) 



naok 2 

£TF(k) + <7o/if>o]. (33) 



Interestingly, the inclusion of correlation, already on 
the linearized LDA level, has introduced the exchange- 
correlation kernel f xc in the corrected TF dielectric func- 
tion, and thus lead to a dependence on local-field correc- 
tions that would ordinarily only be expected in sophis- 
ticated beyond-RPA dielectric functions. 1 ' 2 ' 7 The wave- 
vector dependence of the local-field correction, however, 
is lost in the linearization. 



VI. TF DIELECTRIC FUNCTION IN 
TWO-DIMENSIONAL ELECTRON LIQUIDS 

Next, we extend the above analysis to two dimensions, 
where our starting point is the 2-d local-density approx- 
imation to the kinetic-energy functional, 

T tfM = ^ J d 2 rn{vf. (34) 

(Here and below 2-d energy functionals and dielectric 
functions are distinguished from their 3-d counterparts 
by a subscript 2. For all other quantities no confusion 
can arise, and to avoid cluttering the formulas we use 
the same symbol for them in 2-d as in 3-d.) By follow- 
ing exactly the same steps as before we find the Euler 
equation 



m 



n(r) + v ext (r) + v H (r) - (J, = 0. (35) 



Unlike Eq. (8) in 3-d, this density-potential relation is 
already linear, and no further approximation is needed to 



■5 



write it in the form n(r) = no 



i(r). The 2-d Fermi 



wave vector is related to the 2-d density via n = kp/2n, 
and in terms of the 2-d density of states go — m/irh 2 
one immediately identifies rii n d = —gov s . All remaining 
steps are the same as in 3-d, and from the 2-d Fourier 
transform ^[l/r] — 2n/k one readily obtains 



Utest(k) 



iw(k) 



1 + e 2 g F{l/r} 1 + 



a k 



v s (k) = 

VtestQt) 
£TF,2(k) ' 



For 



screening of a point charge (36) becomes 



V s (k) 



2^e 2 



k + h 



TF 



(36) 



(37) 



where in 2-d k TF = 2/a . 
independent expression 



The resulting density- 



2 1 

e T F,2(k) =1-1 - 

do k 



(38) 



for the 2-d TF dielectric function is well known, e.g., 
in the context of 2-d electron liquids in semiconductor 
heterostructures. 15 ' 18 

Interestingly, the 2-d TF dielectric function follows 
rigorously from the 2-d TF kinetic energy functional, 
without an additional linearization. Hence, in this 
sense, (36) is a stronger result than its 3-d counterpart 
(13). 16 However, this stronger result is an artifact of two- 
dimensionality, arising because the 2-d kinetic energy is 
quadratic in n(r), and its density-derivative linear. 

On the other hand, the fact that erF,2(k) is obtained 
from TxF,2[n\ in the same way eTF(k) is obtained from 
TtfIit] shows that the relation between energy function- 
als and dielectric functions in itself is not an artifact of 
three-dimensionality. 17 



VII. TFD DIELECTRIC FUNCTION IN 
TWO-DIMENSIONAL ELECTRON LIQUIDS 

Finally, we work out the form of the exchange cor- 
rections to the two-dimensional Thomas-Fermi dielectric 
function. The per-particle exchange energy of a uniform 
two-dimensional electron liquid is — 8v / 2/(37rr s )i?y. 19 ' 20 
where r s = l/(a ^/nn) and Ry = e 2 /2ao- The local- 
density approximation for E Xt2 [n] is thus 

E L J A [n] = -\{^e 2 j<?rnV\ (39) 

The 2-d TFD Euler equation can be cast in the same 
form as before, n(r) = n + n ind (r), where in 2-d n = 
kp/2-n- and n ind = -gp^Vs, with v s = v ext +v H + v x and 
v x = —2e 2 y/2n/ii. The Hartree potential arising from 
the unperturbed density uq and the potential energy in 
the positive background (contained in v ext ) cancel again, 



and by means of the convolution theorem the Fourier- 
transformed screened potential becomes 



v s (k) = v ext (k) + e 2 T 2 [n md (r)]T 2 



^ 2 ^-T 2 [^\. 



(40) 

The Fourier transforms in the second term on the rhs 
are easily evaluated by using rii n d = — go, 2 v s (k.) and 
F 2 \\-jr\ = 2n/k. To evaluate the Fourier transforms in 
the third term we linearize in terms of n ind /n a , as before, 
and obtain 



T 2 [Vno + n md \ w T 2 [V^o] 



1 



T 2 [n». 



(41) 



Hence, for the 2-d screening of a point charge with 
vtest(r) = e 2 /r, we obtain 



«.(k) = 



27re 2 



1 



k 



1 ~r a k ir ' 3 



(42) 



where we have used r s = ^/2/aokF for the 2-d density 
parameter r s , and, as before, kS(k) = 0. As in the 3- 
d case, the derivation from a kinetic-energy functional 
allows us to build exchange effects in a simple way into 
the TF dielectric function. 

Interestingly, our final expression 



eTFD,2W = H 1 r s = £TF,2(k) r s 

a k TT TT 



(43) 



has aquired a density dependence, through r s {n) = 
1/aoy/nTT- Inclusion of exchange thus removes the un- 
physical density- independence of £TF.2(k) obtained in the 
2-d TF approximation (36). The behaviour of the 2-d 
TF and TFD dielectric functions as a function of wave 
vector and density is illustrated in Figs. 2 and 3. Cor- 
relation corrections could be included in two dimensions 
in the same way as in three dimensions. In the interest 
of brevity we refrain from showing the straightforward 
extension of the above expressions to this case. 

On the other hand, we stress that in 2-d we obtain 
the same chain of equalities as in 3-d, relating the addi- 
tive corrrection to the dielectric function with multiplica- 
tive renormalization of the sound velocity, compressibil- 
ity and spin susceptibility of the 2-d electron gas: 21,22 



£TF£>,2 (k) £TF,2 (k) = 1~ = 1-^| = 1-^ = - 



X 



V2r s 

TT 

(44) 

Eq. (23) thus carries over to 2-d, with a modified density 
dependence. 



VIII. SUMMARY AND ASSESSMENT 

While the 2d and 3d TF dielectric functions ej , ^(k) 
and eTF,2(k), given in Eqs. (15) and (38), are well known 
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1 1.5 2 2.5 3 3.5 4 4.5 5 



r 

s 

FIG. 3: Density dependence of the 3-d TF (dashed line), 
3-d TFD (full line), 2-d TF (dotted line) and 2-d TFD 
(dash-dotted line) dielectric functions at fixed wave vector 
k = l/3ao (where A = 2n/k — 6nao is the wave length of 
the density variations in n in d(r)). The change of the TF di- 
electric function due to exchange corrections is a downshift 
that is more pronounced at low densities (large r s ). In 2-d 
it is only because of the exchange corrections that a density- 
dependence appears in the dielectric function. 

(although not normally derived from density function- 
als), their TFD counterparts ctfd(^-) and ctfd. 2(k), 
Eqs. (22) and (43), and the correlation-corrected expres- 
sion (33) are apparently derived for the first time in this 
work. Wc do not expect simple models of this type to 
substitute the sophisticated models available from dia- 
gramatic many-body physics, 7 or the material-specific 
dielectric functions that can be obtained from density- 
functional theory. 12 However, each of the four extensions 
of the 3-d TF functional studied here displays interesting 



features that were not obvious from the 3-d TF case, but 
are revealed by deriving model dielectric constants from 
density functionals. 

Specifically in the 3-d TFD case, we found that the 
exchange contribution to the TF dielectric function yields 
a term of the same form as the long- wavelength expansion 
of the RPA, if that expansion is carried out consistently 
in orders of k. The standard way of identifying erF(k) 
as the k — > limit of the static RPA is not consistent in 
orders of k. Further correlation terms lead, already on 
the linearized LDA level, to a simple dependence of the 
dielectric function on local-field factors. 

Unlike the connection between the TF dielectric func- 
tion and the Yukawa potential in real space, that be- 
tween the TF dielectric function and energy functionals 
carries over to two-dimensions, where £TF,2(k) turns out 
to be a more robust expression than eTF(k) in 3-d, be- 
cause it can be derived from the 2-d TF kinetic-energy 
functional without any linearization. Upon including ex- 
change (TFD), a density dependence reappears in the 
2-d dielectric function, which is lost in the usual 2-d TF 
dielectric function. In both 2d and 3d, the additive ex- 
change contribution to the dielectric function turns out 
to be the same as the multiplicative exchange renormal- 
ization factor of other observables, such as the spin sus- 
ceptibility. 

The existence of a connection between a class of energy 
functionals and a set of model dielectric functions opens 
the question whether more detailed models for e(k) can 
be obtained from the more complex energy functionals 
in use in modern density-functional theory. Models for 
frequency-dependent dielectric functions may be accessi- 
ble via time-dependent density- functional theory. 23 
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